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ABSTRACT 

Knowledge of intrinsic shape and orientation of galaxy clusters is crucial to understand their 
formation and evolution. We propose a novel model which uses Bayesian inference to deter- 
mine the intrinsic form of the hot intracluster medium of galaxy clusters. The method exploits 
X-ray spectroscopic and photometric data plus measurements of the Sunyaev-Zel'dovich ef- 
fect (SZe). The gas distribution is modelled with an ellipsoidal parametric profile who can fit 
observed X-ray surface-brightness and temperature. Comparison with the SZ amplitude fixes 
the elongation along the line of sight. Finally, Bayesian inference allows us to deproject the 
measured elongation and the projected ellipticity and constrain the intrinsic shape and orien- 
tation of the cluster. We apply the method to the rich cluster Abell 1689, which was targeted 
by the Chandra and XMM satellites as well as by several SZe observatories. Observations 
cover in detail a region < 1 Mpc. Our analysis favours a mildly triaxial cluster with a mi- 
nor to major axis ratio of 0.70 ± 0.15, preferentially elongated along the line of sight, as 
expected for massive lensing clusters. The triaxial structure together with the orientation bias 
can reconcile X-ray with lensing analyses and supports the view of A1689 as a just slightly 
over-concentrated massive cluster not so far from hydrostatic equilibrium. 
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1 INTRODUCTION 

Clusters of g alaxies are the most recent bound structures to form in 
the Universe ( Voit 2005). Their intrinsic shape contains evidence of 
the assembly process of structures. Accurate knowledge of forms 
is important on its own by probing the cosmic structure formation. 
It suggests how material a ggregates from large-scale perturbations 
dWestl(l994 IJing & Sutd l2002) and contains evidence about the 
nature and mechanisms of interaction of b aryons and dark matter 
dLee & Sutdl2003l:lKazantzidis et alj2004h . 

The complex structure of halos also affects the estimation of 
quantities which are crucial in any attempt at high precision cos- 
mology. Favouring a simple spherical model over more realistic 
shapes might c ause a significant bias in estimating the cluster mass 
dGavazzill2005h. the in ner matter density slope and the concentra- 
tion dOguri et alj|2005ri . 

Accurate knowledge of the intrinsic structure is also criti- 
cal when comparing observations with theoretical predictions. The 
observed concentration-mass relation for galaxy clusters has a 
slope consistent with what found in TV-body numerical simulations, 
though the normalization factor is higher dComerford & Nataraiad 
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1 20071 : [Ettori et al] |2010l) . Disagreement between theory and obser- 
vation might be explained by orientation and shape biases. In fact, 
triaxial halos can be much more efficient lense s than their more 
spherical counterparts (Oguri&Blandford 2009) with the strongest 
lenses in the Universe expected to be a highly biased population 
preferentially orientated along the line of sight. 

Classical attempts to determine intrinsic three dimensional 
forms were based on statistical approaches consi sting in t he in- 
version of the d i stribution of apparent shape s (Hubble: 1 192 



Binggeli 1980; Binnev & de Vaucouleursl 



Fasano & Viol 199ll; Ide Theiie et alj Il995l; iMohr et alj 



1991 : 
liooil: 
ct al. 



Basilakos et al. 2000; Cooray 2000; Thakur & Chakraborty 
Alam & Rvdenll2002l ; lR"vdedll996l ; Wonis et alj|2004l ; IPaz 
l2006t) . With the exception of disc galaxies, either prolate-like or 
triaxia l shapes appear to dominate all cosmic structure on a large 
scale. Kawahara (2010) found that the observed probability density 
function of the projected axis ratio of a sample of X-ray clusters 
was compatible with that of a population of triaxial halos in hydro- 
static equilibrium consistent with TV-body simulations. 

The d etermination of the f orm of single ob j ects is a more re- 
cent topic dZaroubi et alii 19981 ; lReblinsk\ll200d ; iDore et aljfeooil : 
IPuchwein & Bartelmandl2006l) . Clusters of galaxy are very inter- 
esting targets since they can be probed with very heterogeneous 
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Figure 1. Exposure-corrected image of one of the Chandra observation 
used for the spatial analysis. The full lines are the best-fit ellipses enclosing 
50, 60, 70 and 80 per cent of the light for a fixed X-ray centroid. 



data-sets at very different wave-lengths from X-ray surface bright- 
ness and spectral observations of the intra-cluster medium (ICM), 
to gravitational lensing (GL) observations of the total mass distri- 
bution to the Sunyaev-Zel'dovich effect (SZe) in the radio-band. 
On the theoretical side, it has be e n ascertained that the de projec- 
tion is not unique tevbickjH987l : iGerhard & Binnevlll996h . Even 
assuming the cluster to be a triaxial ellipsoid, the only quantity that 
can b e univocally determined is the elongation along the line of 
sig ht dSerenoi r2007). This allows to break the degeneracy with the 
distance and to g et an unbiased estimate of the Hubble constant 
dFox & Penll2002h . 

On the observational side, only a few works have tried to in- 
fer shape or orientation of single objects. Combined use of X-ray 
and SZe data allows to constrain the shape of the ICM without any 
assum ption regarding eq u ilibriu m or geometry. |Pe Filippis et alj 
J2005I) and ISereno et alj (2006) performed a parametric analy- 
sis of a sample of 25 clusters, finding that prolate rather than 
oblate shapes seem to be preferred, with signs of a more general 
triaxial morphology. Mah davi & Chand feOllh derived a model- 
independent expression for the minimum line-of-sight extent of the 
hot plasma and applied it to the Bullet Cluster. A different approach 
is based on lensing observations. Surface density maps of the total 
matter distribution from weak or strong lensing have been depro- 
jected explo i ting some a priori assumptio ns on their intrin sic shape s 



Qguri et alj d2005h; ICorless et al.1 j2009h : ISereno et alj d2010al 
Sereno & Umetsul ( 1201 lh . 

Here, we propose a parametric method to infer the intrinsic 
shape and orientation of clusters based on deep X-ray and SZe ob- 
servations. First, we model the gas density and fit to the data to 
obtain a direct measurement of the elongation along the line of 
sight. Then, we exploit Bayesian methods which enable to make 
an inference about a number of variables larger than the number of 
the observed data. This allows to constrain both intrinsic shape and 
orientation. 

The method is applied to Abell 1689 (A1 689), see Fig. [T] 
a ver y luminous cluster at r edshift z = 0.183 (Bro adhurst et al.1 
120051 : I Limousin et "al] |2007r) . Strong lensing analyses have pro- 
vided a consistent picture of the mass distribution in the in- 



ner < 300 kpc regions dBroadhurst et al] 120051 ; lHalkola et all 
200 6; Limousin et al j2007lJcoe et alj2010h favouring a quite con- 
centrated mass distribution. On the other hand, some disagree- 
ment persists on the larger virial scale. Different weak lens- 
ing analyses suggest somewh a t different degrees of concentration 
dUmetsu & B roadhurs ll |2008| ; lUmetsu et all 120091 ; ICorless et ail 
2009; Limousin et alj |2007). There is some conflict between X-ray 
and lensing analyses, with lensing masses exceeding estimates de- 
rived under the hypothesi s of hydrostat i c equ ilibrium by 30-40 per 
cent in the inner regions. iLemze et alj d2009h combined Chandra 
X-ray brightness measurements and joint strong/weak lensing mea- 
surements under the hypothesis of hydrostatic equilibrium. They 
found that the resulting equilibrium temperature exceeds the ob- 
served temperature by 30 per cent at all radii. 

Disagreement can be reduced properly accounting for devi- 
ations from either spherical symmetry or hydrostatic equilibrium. 
The level of hydrostatic equil i brium in A 1689 has been recently 
reassessed. iKawaharada et al.1 d2010h found regions with low gas 
temperatures and en tropies deviating from hydrostatic equilibrium. 
Moln ar et a showed that in the core region of relaxed clus- 

ters alike A 1689 a significant non- thermal pressure support of ~ 20 
per cent is originated by subsonic random gas motions within 1-10 
per cent of the virial radius. 

A number of triaxial lens i ng analyses of A 1689 has 
been performed dOguri et alj|200l ; ICorless et ai]|2009h . Recently, 
ISereno & Umetsu < 201 1 ) developed a method for a full three- 
dimensional analysis of strong and weak lensing data. They found 
evidence for a mildly triaxial lens (minor to major axis ratio ~ 
0.5±0.2) with the major axis orientated along the line of sight. The 
triaxial shape is compatible with a halo slightly over-concentrated 
but still consistent with theoretical predictions. IPeng et alj ( 2009) 
suggested that a prolate ellipsoidal configuration for the gas distri- 
bution, aligned with the line of sight and with an axis ratio of ~ 0.6, 
could solve the central mass discrepancy betw een lensing and X- 
ray estimates. Recently, Ifvlorandi et alj d20 1 lh combined lensing 
and X-ray data under some very restrictive hypotheses. They as- 
sumed the cluster to be aligned with the line of sight and shaped 
deviations from hydrostatic equilibrium with a very peculiar mod- 
elling. They also exploited results from numerical simulations re- 
lating the shapes of gas and matter distribution, which are affecte d 
by quite large errors for very elongated clusters dLee & Sutol2"003l) . 
They found an axial ratio for the matter distribution of ~ 0.5. 

Here, we address the shape shape determination of A 1689 
with a novel model. The paper is organised as follows. In Sec. [2] 
we discuss the relation between intrinsic and projected quantities of 
triaxial ellipsoids. Section[3]hsts some measurable X-ray and SZe 
quantities. Observations and data reduction are described in Sec. [4] 
In Sec. [5] we detail how the gas distribution was modelled and fit- 
ted to the data. Section[6]lists some sources of systematic errors. In 
Sec.|7]we perform the deprojection. Section[8]is devoted to some fi- 
nal considerations. Throughout the paper, we assume a flat ACDM 
cosmology with density parameters Q.m = 0.275, = 0.725 and 
Hubble constant Hp = WOh km s -1 Mpc - \ h = 0.702 ± 0.014 
dKomatsu et alj|201 il) . At the A1689 distance, 1" corresponds to 
3.08 kpc. 



2 TRIAXIAL ELLIPSOIDS 

High resolution A^-body simulations have sh own that the den- 
sity profiles of m assive halos are aspherical dJing&Sutoll2002l; 
iLee & Sutol |2003|) . The electron density n of the intra-cluster 
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Figure 2. The orientation of a triaxial ellipsoid. The z^ nt -coordinate axes 
are oriented along the principal axes of the ellipsoid and define the intrinsic 
reference system. The x° bs -axes defines the observer's frame. The line of 
sight is oriented along the Xj bs -axis. The Euler's angle i9, ip and i/j are 
displayed. Thick lines denote the ellipsoid axes. l s is the major axis. 





if" 

Figure 3. The ellipsoid of Fig. [2] as seen from the observer. Left panel: 
projection of the ellipsoid into the plane of the sky. The observer sees an 
ellipse with major axis l p . The axes of the ellipse are plotted as thick lines. 
The dashed line is the projection in the plane of the sky of the Xg nt -axis. 
Right panel: ellipsoid as seen from above. This would be the view of an 
observer perpendicular to the line of sight, i.e., located on the xj bs -axis. 
'l.o.s. i s me na h° size along the line of sight (z2 bs -axis). 



medium (ICM) can be assumed to be constant on a family of simi- 
lar, concentric, coaxia l ellipsoids. Thes e assumptions were also ob- 
servationally verified ( Kawahara 2010). 

The gas distribution in clusters of galaxies in hydrostatic equi- 
librium traces the gravitational potential. Since we are considering 
a triaxial elliptical gas distribution, the gravitational potential turns 
out to be ellipsoidal too. These gravitational potential can turn un- 
physical for extreme axial ratios, giving negative density regions or 
very unlikely configurations, but as far as either inner regions or 
small eccentricities are considered, they can provide very suitable 
approximations. 

An ellipsoidal ICM profile can be expressed as a function of 
only one radial variable £, which runs along the major axis and 
replaces the spherical radius. The minor (intermediate) to major 
axial ratio is denoted as qi (q 2 ) with < qi ^ q% ^ 1; we also 



use the inverse ratios, < = l/q% ^ 1. For a prolate shape, 
qi — 92 ^ 1; a oblate ellipsoid has qi ^ q% = 1. 

Three Euler's angles, i9, p and ip, relate the intrinsic to the 
observer's coordinate system, see Fig. [2] The angle -d quantifies the 
inclination of the major axis with respect to the line of sight. 1? 
and ip fix the orientation of the line of sight in the intrinsic system. 
The third angle ip determines the orientation of the cluster in the 
plane of the sky. Since we are interested in the cluster shape, we 
can neglect this angle in the following. 

When viewed from an arbitrary direction, quantities co nstant 
on similar ellipsoid s project themselves on similar ellipses dStarkl 
1 19771 : ISeren3l2007h . see Fig. [3] The ellipticity and the orientation 
of the projected ellipses depend only on the intrinsic geometry and 
orientation of the system. The axial ratio of the major to the minor 
axi s of the observ ed projected isophotes, e p (js 1), can be written 
as teinggelilll980h . 



3+1+ V(j - I) 2 + 4fc 2 



(1) 



j + l- V(J - 2 + 4fc2 
where j, k and I are defined as 

• — 2. - 2- 1 * 2 2n 2 |2 2 n ■ 2 /^\ 

3 — e i e 2 sm v + e 1 cos v cos <p + e 2 cos v sin p>, (2) 



k — (e 1 — e 2 ) sin ip cos <p cos 1}, 
I = e 2 sin 2 tp + e 2 cos 2 <p. 



(3) 
(4) 



We also use the ellipticity e = 1 — 1 /ep. 

Ellipsoids map into ellipses. Let l s be the intrinsic major axis 
of the ellipsoid, see Fig. [2] The corresponding projected ellipse on 
the plane of the sky has a major axis l p , see Fig. [3] 

1/2 







eie 2 



.f 



1/4 



(5) 



where / is a function of the cluster shape and orientation, 

2 2 2 2 2 2 2 

f — ex sin 1? sin <p + e 2 sin v cos ip + cos 

The quantity 
h.o.s. = Is/ v/ 

is the half-size of the ellipsoid along the line of sight, i.e., as 
seen fro m above, see Fig. [3] Projection and elongation are related 
through dSerendl2007h 



(6) 



(7) 



^l.o.s. — 



h_ 

eA ' 



(8) 



which represents the definition of the elongation eA- In terms of 
intrinsic parameters, 

'a = 1 I 17 '/ 374 . (9) 



ep 
eie 2 



eA is the ratio between the major axis of the projected ellipse in the 
plane of the sky and the size of the ellipsoid along the line of sight. 
It quantifies the elongation of the triaxial ellipsoid along the line of 
sight. If eA < 1, then the cluster is more elongated along the line 
of sight than wide in the plane of the sky, i.e., the smaller eA, the 
larger the elongation along the line of sight. 



3 PROJECTION 

The projected map F20 of a volume density F30, which is constant 
on surface s of constant ellipsodial radius C, is elliptical o n the plane 
of the sky dStarkll977l : rSerendl2007l : [Sereno et al|2010al) . F 2 u has 
the same functional form of a spherically symmetric halo, 
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e A J Xf 



Pi) 



(10) 



where = f/Z 8 is the integration variable, £ is the elliptical ra- 
dius in the plane of the sky (measured along the projected major 
axis), = £/Zp, l s is the length scale of the 3D density, Zp is its 
projection on the plane of the sky, and pi are the other parameters 
describing the intrinsic density profile (inner or outer slope, tidal 
radius, ...)• If we know the form of the intrinsic profile F3D, the 
only unknown in the projection is the undetermined geometrical 
factor e7 . 



3.1 X-ray Surface Brightness 

Cluster X-ray emission is due to bremsstrahlung and line radia- 
tion resulting from electron-ion collisions in the high temperature 
plasma (fcp,2o ~ 8-10 keV, with fcs being the Boltzmann con- 
stant). The observed X-ray surface brightness (SB) can be written 

as 



SB° bs = 



n 2 e A* a (T e ,Z)dl, 



(11) 



4*(l + z)» JLOmt 

where T c is the intrinsic temperature, Z the metllicity and Af 5 is 
the effective cooling function of the ICM in the cluster rest frame, 
which depe nds on the ener gy-dependent effective area of the in- 
strument dReese et alj|2010h . We comp uted predicted emissivities 
through the XSPEC software package dAmaud[l996l) adopting an 
optically-thin plasma (mekal) emission model . 



3.2 The Sunyaev-Zel'dovich effect 

Photons of the cosmic microwave background (CMB) that pass 
through the hot ICM of a cluster interact with its energetic 
electrons through inverse Compton scattering, slightly distort- 
ing th e CMB spectrum. This is t he Sunyaev-Zeldo vich effect 
(SZe) iSunvaev & Zeldovichl Il970l : iBirkinshawl 1 19991) , which is 
proportional to the electron pressure integrated along the line of 
sight. The amplitude of the signal is proportional to the Compton-y 
parameter, 



y 



CTT^E 



nTdL 



(12) 



where <tt is the Thompson cross section, m c the electron mass, c 
the speed of light in vacuum. 

The measured temperature decrement ATsz of the CMB for 
an isothermal plasma is given by 



AT SZ = /sz(i / ,T)Tcmb2/ 



(13) 



where Tomb is the temperature of the CMB and fsz{v,T) ac- 
counts for relativistic corrections at frequency v. 

An overall measure of the thermal energy content in a cluster 
is given by Y, which is the Compton- y parameter integrated from 
its center out to a projected radius ta, 



Y A = 



ydQ, 



(14) 



where JIa is the solid angle of the integrated patch covered within 
a radius r a. Y is a more robust quantity for observational tests than 
the central value of Compton-y, because it is l ess dependent on th e 
model of gas distribution used for the analysis teenson et alj2004T) . 



Table 1. Measured ellipticities from Chandra data. The rotation angle 6 C 
is measured counterclockwise. The Position Angle PA (= 6 e — 90 deg) is 
defined equal to zero at North. 



Exposure ID 


e 


0e 


6930 


0.152 ±0.036 


-77.7 ±5.3 


7289 


0.148 ±0.038 


-78.3 ±3.6 



In addition, integrating the Compton-y out to a large projected ra- 
dius diminishes (though does not completely remove) effects re- 
sulting from the pr esence of strong ent ropy features in the central 
regions of clusters (Mc Carthy et al.ll2003l) . 



3.3 Observed temperature 

The temperature T c in Eqs. QT] and [12] is meant to be the intrinsic 
temperature. The spectroscopic temper ature Tp measured by space 
observatories is well approximated by (Ma zzotta et al. , 2004) 



T P = 



/ WT c dV 
JWdV ' 



W = 



3/4- 



(15) 



4 DATA 

4.1 X-ray surface brightness 

We have reduced the Chandra data, see Fig. Q] using the CIAO 
data analysis package - versi on 4.3 - and the asso ciated calibra- 
tion database CALDB 4.4.1 jFruscione et"ai1l2006h . Out of the 5 
exposures available, we have considered the 2 longest ones, Obs 
ID 6930 (nominal expsoure time of 80 ks) and 7289 (80 ks), both 
in ACIS-I configuration. The level- 1 event files were reprocessed to 
apply the appropriate gain maps and calibration products and to re- 
duce the ACIS quiescent background. The standard event selection 
has been applied through the CIAO routine acis_process_events in 
"very faint" mode configuration. After the exclusions of extended 
and point sources present in the field and selected with the vtpde- 
tect routine, a cleaning of the background light curve has been per- 
formed with the deflare routine for a total exposure good time of 
75.6 ks (Obs ID 6930) and 75.1 (Obs ID 7289). 

The images in the 0.7-2 keV band have been corrected for the 
exposure map, and then for the presence of the point-sources using 
the CIAO routine dmfilth to produce the image of only extended 
emission associated to the cluster. A centroid is defined within the 
source detected with vtpdetect and related to the cluster ICM using 
the CIAO routine dmstat. The centroid is then fixed and the routine 
dmellipse is run over the exposure-corrected image of the cluster. 

Measured ellipticity e and orientation angle are reported in Ta- 
ble!]] The best-fit results associated to the fit enclosing 80 per cent 
of the cluster light are considered, see Fig. [T] Being the two esti- 
mates consistent, we fix ellipticity and angle to the measurements 
obtained in the exposure 7289. Our estimate of th e ellipticity iscon- 
sistent with the value of 0.16 ± 0.04 measured in Kawahara (2010) 
exploiting data from XMM. 

The surface brightness profiles are then extracted in elliptical 
annuli. Finally, the two profiles are combined in a single profile by 
summing the counts and propagating the error on the estimate of 
the background. We evaluate the extension of the detectable ICM 
emission by requiring that the signal-to-noise ratio is larger than 2. 
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We measure a maximum radius of 5.9' (1093 kpc). The resulting 
surface brightness profile is plotted in Fig.[4] 

4.2 X-ray spectroscopic analysis 

The Observation Data Files (ODF) from the XMM-Newton obser- 
vation have been processed to produce calibrated event files using 
the XMM-Newton Science Analysis System (SAS vlO). 

To search for periods of high background flaring, light curves 
for pattern ^ 4 in the 10-12 keV have been produced for two MOS 
detectors. The soft proton cleaning w as performed using a doubl e 
filtering process similar to the one o f ILeccardi & Mole ndi (2008). 
We extracted a light curve in 100s bins in the 10-12 keV energy 
band by excluding the central CCD, applied a threshold of 0.20 cts 
s _1 , produced a GTI file and generated the filtered event file ac- 
cordingly. We then extracted a light curve in the 2-5 keV band. 
The residual intervals of very high background were removed us- 
ing a 3a clipping algorithm and the light-curves were then visually 
inspected to remove the background flaring periods not detected by 
the algorithm. The clean exposure times for the two MOS detectors 
were ~ 36 ks. PN data were not used in the analysis because of 
problems in the background modelling and inconsistencies in the 
determination of the temperatures with the MOS detectors. 

The background in our MOS spectra has been modelled 
(instead of subtracted) fo llowing the procedure developed by 
ILeccar di & Molendi (2008). This procedure is preferable respect 
to a direct background subtraction because it allows us to use Cash 
statistics and avoid problems due to the vignetting of the back- 
ground spectra which would need to be extracted at larger off-axis 
angle than the cluster spectra where the vignetting of XMM mir- 
rors is not negligible. The main aspects of this procedure are de- 
scribed below, further details can be found on their paper. The back- 
ground parameter were first estimated in a region free of cluster 
emission (10'-12' annulus). The background model considered in- 
cluded the thermal emission from the Galaxy Halo (HALO, XSPEC 
model: apec), the cosmic X-ray background (CXB, XSPEC model: 
pegpwrlw), a residual from the filtering of quiescent soft protons 
(QSP, XSPEC model: bknpower), the cosmic ray induced contin- 
uum (NXB, XSPEC model: bknpower) and the fluorescence emis- 
sion lines (XSPEC model: gaussian). The Response Matrix File 
(RMF) of the detector was convolved with all the background com- 
ponents, while the Ancillary Response File (ARF) was convolved 
only with the first two components (HALO and CXB). The normal- 
ization of the QSP component was fixed at the value determined 
from measuring the surface brightness SB ln in the 10'-12' annu- 
lus, and comparing it to the surface brightness SB out calculated 
outside the field of view in the 6-12 keV energy band. Since soft 
protons are channeled by the telescope mirrors inside the field of 
view and the cosmic ray induced background covers the whole de- 
tector, the ratio Rsb = SB- m /SB on t is a good indicator of the 
intensity of residual soft protons and was used for the modelling 
(Norm(QSP) = 0.03(1 - Rsb))- We determined also the la 
error for the background parameters to be used in the fit of the 
cluster spectra. To fit the cluster spectra, the background param- 
eters (and their la errors) were then rescaled by the area where 
the cluster spectra were ex tracted. Appropriate correction factors 
(Leccardi & Molendi 2008, dependent on the off-axis angle) were 
considered for the HALO and CXB components and a vignetting 
factor (corresponding to 1.858 — 0.078r, with r the distance from 
the center of the annulus) was applied to the QSP component. All 
these rescaled values (and their la errors) were put in a XSPEC 
model having the same background components used in the fit of 



Table 2. Integrated Compton parameter within 600 kpc, in units of 
10 _10 rad 2 , as estimated by different observatories. uii ns is a statistical 
weight. References are listed in column 4. 



Observatory 


^2500 


Wins 


reference 


BIMA/OVRO 


2.0 ±0.4 


1/2 


Bonamente et al. (2006) 


BIMA/OVRO 


2.8 ±0.5 


1/2 


Reese et al. (2002) 


AMIBA 


2.8 ± 1.1 


1 


Huang et al. (20101 


SuZIE 


3.6 ± 1.2 


1/2 


Holzapfel et al. (19971 


SuZIE 


4.4 ± 1.2 


1/2 


Benson et al. (20041 


WMAP 


2.3 ±0.7 


1 


Lieu et al. (20061 


SKA 


2.3 ±0.8 


1 


Grallaet al. (2011) 


SCUBA 


5.8 ±2.4 





Zemcov et al. (2007) 



the 10'- 12' annulus plus a thermal mekal model for the emission 
of the cluster having the temperature, the abundance and the nor- 
malization free to vary. The la errors on the parameters were used 
to fix a range where the normalizations of the background com- 
ponents are allowed to vary. This model was used in addition to 
the one temperature thermal model used to fit all the spectra and 
described below. 

We divided the emission in five concentric elliptical annuli. 
The temperature profile is plotted in Fig. [6] Each annulus spectrum 
was analyzed with XSPEC vl2.5.1 lArnaudl 19961) and fitted with a 
single-temperature mekal model ( ILiedahl et all 19951) with Galac- 
tic absorption (wabs model), in whic h the ratio between t he ele - 
ments was fixed to the solar value as in lAnders & Grevessel (l989). 
The fits were performed over the energy range 0.7 — 8.0 keV for 
both of the two MOS detectors. We fitted simultaneously MOS1 
and MOS2 spectra to increase the statistics. The free parameters 
in our spectral fits are temperature, abundance, and normalization. 
Local absorption is fixed to the Gala ctic neutral hydrogen column 
density, as obtained from radio data dKalberla et al]|2 005), and the 
redshift to the value measured from optical spectroscopy. We used 
Cash statistics applied to the source plus backgrounqj, which re- 
quires a minimal grouping of the spectra (at least one count per 
spectral bin). 

4.3 SZe 

We reviewed available SZe analyses of A 1689 in literature, see Ta- 
ble^ From published profiles and amplitudes, we estimated the 
parameter Y in a circular region of outer radius 600 kpc. The 
mean over-density within this ra dius with respect to th e cosmic 
background density is ~ 2500 dBonamente et al J 1 20061) . Due to 
angular resolution and instrument sensitivity, analyses exploiting 
either BIMA/OVRO, AMIBA, SuZIE, WMAP and SCUBA data 
usually refer to such radius to determine the SZe amplitude. On the 
other hand, the in tegrated Compton parameter from the SKA data 
(Gralla et alj|201ll) has been provided within ~ 500", i.e., within 
an overdensity radius of 500, way larger than the X-ray spectro- 
scopic coverage of the cluster. In order to have a coherent combined 
X-ray plus SZe analyse s, we rescaled the i r resu lt within 600 kpc 
assuming the profile in iBonamente et al] ([2006). Relativistic cor- 
rections were applied when necessary. 

Since each analysis assumed a fixed parameterization for the 
density and temperature profile, usually the isothermal /3-model, 
we added an error of 10 per cent accounting for dependence on 

1 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XSappendixCash.html 
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Table 3. Contributions to the total \ 2 value under different hypothesis. We 
consider: ij = and <5SB syst = (col. 3); r) > and 5SS syst = (col. 
4); r/ = and 6SB s y st = 10" 3 (col. 5). 





N data 


ssB s y st = o 
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ass 8 ^ = lo- 3 


Xx+sz 


68 


422.7 


422.7 


80.7 


XSB 


56 


415.9 


415.9 


73.9 
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4.3 


4.3 


4.2 


x\ 


7 


2.5 


2.5 


2.5 



modelling (Benson et al. 2 0041 ; iLiao et al.l l2010). This systematic 
error, which was added in quadrature, was usually bested by the sta- 
tistical uncertainty inferred with standard error propagation, so that 
its exact value does not matter significantly. Most of the analyses 
exploited fixed scale radii and slopes for the density profiles from 
independent X-ray analyses. This approach underestimates param- 
eter degeneracy and statistical uncertainties, so we added in quadra- 
ture a further error of 2 5 per cent. The exception is the analysis in 
Bonamente et al. I JUS) which fitted the scale radius together with 
the amplitude to the visibility map. The additional error in this case 
is not needed. 

Some data sets were analysed by independent groups. We 
considered all published results and weighted each analysis in 
such a way that the statistical weight of each instrument is the 
same, see the u>i ns t parameter in Table [2] Finally, measured SZ 
amplitudes from SCUBA are consistentl y higher than resu lts in- 
ferred from low-frequency measurements JZemcov et alj2007h . We 
did not consider the SCUBA result in our investigation. The fi- 
nal weighted mean for the integrated Compton parameter within 
600 kpc is I2500 = 2.5 ± 0.4. 
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Figure 4. Radial profile along the projected major axis of the 0.7-2 keV 
surface brightness measured by Chandra (points with error bars) in elliptical 
annuli. The full line plots the surface brightness profile predicted by the best 
fit model. 



5 THE GAS DISTRIBUTION 



£ [kpc] 



We model the 3D electronic density and temperature with paramet- 
ric profiles. Distributions are assumed to be coaligned and ellip- 
soidal, with constant e ccentr i city and orientation. Profiles are taken 
from IVikhlinin et all d2006h : lEttori et alj J2009l) . For the density 
profile in the intrinsic coordinate system, we use 



n = n | — 







V/2-3P/2 


H0 


2 
3 



















,(16) 



where no is the central electron density, r c is the core radius, r t (> 
r c ) is the tidal radius, fi is the slope in the intermediate regions, r\ 
is the inner slope and (0 ^)7(^ 2.5) determines the outer slope. 
For the temperature profile, we use 



T = T 



(CAt) 



(CAcooir- 01 + Tmin/Tb 



[1 + (C/r T ) f, T] CT/i " r (CAcool) a ™°' + 1 



(17) 



where the radius r coo i and the slope a CO oi characterise the central 
cool core with central temperature T m i n ; To is the temperature in 
the intermediate region; the radius tt and the slope ct describe 
a decrement at large radii, with the width of the transition region 
fixed by &t. 

The metallicity is fixed to the mean observed value. 

We are now in position to compare observations with theoret- 
ical predictions. Values of X-ray surface brightness, spectroscopic- 
like temperature and SZ decrement for a given set of parameters 
can be calculated by first plugging the expressions for n, Eq. J 1 6b . 
and T, Eq. dl7K in the corresponding line of sight integrals, and 
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Figure 5. Radial profile along the projected major axis of the residuals 
ASB between observed and predicted surface brightnesses. Each scatter 
is in units of the corresponding observational error SSB. 



then by evaluating the intergrals according to the prescription in 
Eq. [To] Observations and predictions can then be compared with a 
X 2 function, 



Xx+sz 



2 .2.2 

= XSB +XT+XY 



(18) 
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Electronic density profile Temperature profile 

e A n.o[xlO~ 4 cm _3 ] r P [kpc] r t /r c @ 7 T [KeV] r T /r c 

0.66(.53) ± 0.21 0.31(.28)±0.05 47.2(47.1) ± 1.0 13.8(13.9) ± 1.5 0.464(.464) ± 0.005 1.7(1.7) ±0.4 9.79(9.82) ± 0.19 16.4(16.5) ± 1.6 

Table 4. Inferred model parameters for the electronic density profile and the temperature profile. 

. Central location and dispersion are computed as mean and variance of the PDF. Maximum likelihood values are reported in parentheses. 
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Figure 6. Radial profile along the projected major axis of the projected 
temperature measured by XMM (points with error bars) in elliptical annuli. 
Diamonds denote the spectroscopic-like temperatures predicted by the best 
fit model. The projected spectroscopic-like temperature is a function of the 
elliptical radius £ . The full line represents the temperature of the gas distri- 
bution T c as measured along the ellipsoidal radius f , 
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Figure 7. Posterior probability distribution of the elongation e a . 
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Since we are fitting directly to the projected quantities, instead of 
the core radius r c we determine the projected radius rp. They are 
related through Eq. (3)- The remaining scale lengths are determined 



in units of r c . The values of Nsb, Nt and Nsz are reported in 
Tabled 

The models in Eqs. J 16b and J17t have great functional free- 
dom and can adequately describe almost any type of smooth 
profile for extended radial ranges. On the other hand, over- 
parameterization could bring modelling of unrealistic features. The 
projected temperature was measured in 5 elliptical annuli, whereas 
the profile in Eq. dl7t has 8 free parameters. Firstly, since we did 
not notice a decline in the inner region, we fixed T m i n = To. Sec- 
ondly, degeneracy issues do not allow us to disentangle the effects 
of either slope ct or tidal radius rr. We then fixed ct = 0.9. 
Thirdly, the extent of the transition region can not be constrained 
and we used 6t = 2. We end up with only 2 free parameters for 
the temperature profile (To and Tt). 

In order to find which parameters in the modelling of the den- 
sity profile were really needed to reproduce the data, we investi- 
gated the Xx+sz under different assumptions and applie d some 
simpl e information criteria, such as either BIC or AIC (Liddlc 
2007). Minima of the Xx+sz were found through a standard down- 
hill simplex method. Since the value of the minimum did not 
change for either a free r\ or the case 77 = 0, see Tablef3] we found 
no need for an inner power-law-type cusp and set 77 = 0. We are 
then left with 5 free parameters for the density profile. The final 
set of parameters to be determined is then no, rp, /3, r t /r c and 7 
for the density profile, To and vt /r c for the projected temperature, 
plus an overall elongation parameter eA. 

Statistical errors on the surface brightness measurements were 
likely underestimated, see Table f5] At the maximum likelihood 
value, we got Xx+sz — 423, with Xsb — 416 (for 56 points). To 
check if this affects the results, we considered the effect of a larger 
error on the surface brightness by adding in quadrature a systematic 
error of 10~ 3 to each measurement in the SB profile. As a conse- 
quence, at the maximum likelihood the x%b{— 74) got a more 
realistic value whereas the values of xt an d Xy at me minimum 
were not affected. The best-fit set of parameters was negligibly af- 
fected (with a shift < 0.2 per cent for eA and variations < 2 per 
cent on other parameters). We can than conclude that a larger error 
on the SB does not affect our analysis. In fact, best fit values does 
not change significantly whereas the final estimated error on the 
elongation, the parameter we are more interested in, is dominated 
by the statistical uncertainty on the SZe amplitude. 

The best fit model, whose model parameters are summarized 
in the bracketed values in Table [4] reproduce well the data. Ob- 
served values of the surface brightness are plotted in Fig. [4] versus 
the predictions of the best fit model. The scatter plot is in Fig. [5j 
where the gap ASB between model prediction and observed value 
at each radius is represented in units of the observational error 
8SB. Most of the values derived with the best fit model are within 
la from the corresponding measured value, and nearly all of them 
are within 2a. 

Projected temperature profiles are plotted in Fig. [6] The mea- 
sured spectroscopic temperatures are very well reproduced by the 
best fit model. We also plot the gas temperature T . Whereas the 
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spectroscopic-like temperature is plotted as a function of the pro- 
jected elliptical radius £, the intrinsic temperature T a refers to the 
three-dimensional gas distribution and is a function of the ellip- 
soidal radius f . 

The integrated Compton parameter within 600 kpc predicted 
by the model is Y — 2.5, which reproduces perfectly the weighted 
mean of the observed data. This is expected in our triaxial approach 
since the free parameter eA is responsible for the good fit of X-ray 
and SZ data at the same time. 

To asses realistic probability distributions for the parameter 
we performed a statistical Bayesian analysis. The Bayes theorem 
states that 



p(P|d) cx £(P|d)p(P), 



(20) 



where p(Pjd) is the posterior probability of the parameters P 
given the data d, £(P|d) is the likelihood of the data given the 
model parameters and p(P) is the prior probability distribution for 
the model parameters. 

As likelihood, we used Cx+sz oc exp {— Xx+sz/2} ■ As 
priors, we used fiat probability distributions. We explored the 
parameter space by running four Markov chains. We checked 
for chain convergence by verifying that the standard var(chain 
mean)/mean(chain var) indicator was less than 1.2. Results are 
summarised in Table [4] Central location and dispersion for each 
parameter are the mean and the variance, respectively, of the corre- 
sponding marginalised posterior probability function (PDF). Reas- 
suringly, estimates based on the maximum likelihood investigation 
are fully compatible with the results from the Bayesian analysis. 

The posterior probability distribution for the elongation is 
plotted in Fig. [7] The distribution is peaked in correspondence of 
the maximum likelihood value with a long tail at large values, cor- 
responding to rounder configurations. The mean of the distribution, 
see Table[4] is then larger than the best-fit value, even if still com- 
patible within la confidence level. 



6 ADDITIONAL SOURCES OF ERROR 

Other sources of statistical and systematic uncertainty that affect 
the measurement of the elongation have to be considered together 
with the photon-counting statistical uncertainties of the X-ray im- 
ages and spectra, and the statistical uncertainty of the SZe obser- 
vations, which determined the statistical error on eA estimated in 
the previous section. The elongation enters in the equations as a 
overall factor of proportionality. Additional sources of error are 
then similar to what already d iscussed in detail f or the determina- 
tion of the Hubble parameter dReese et alj [2002; Bonamen te et al.l 
l2006j : lDe Filippis et al.ll200"^) . Among the statistical contributions, 
important roles are played by: the uncertainty in Galactic Nn (~ 1 
per cent), the small-scale clumps in the intra-cluster gas (> per 
cent), SZe point sources in the field (~ ±8 per cent), kinetic SZ 
effect (~ ±8 per cent) and CMB anisotropies (< ±2 per cent), 
contribution from the X-ray background (~ ±2 per cent). Finally, 
recent analyses estimated the error on the Hubble constant to be 
~ ±2 per cent JKomatsu et alj|201 ll) . which translates in a similar 
error on eA ■ 

Systematic contributions come from the presence of radio ha- 
los and relics (~ ±3 per cent), the X-ray absolute flux calibration 
(~ ±5 per cent) and X-ray temperature calibration (~ ±7.5 per 
cent). The absolute calibration of a single instrument for SZe detec- 
tion brings about an uncertainty of < ±8 per cent. Since we con- 
sidered several observatories, the corresponding error is reduced to 



91 


92 


COS 1? 


0.70 ±0.15 


0.81 ±0.16 


0.70 ± 0.29 



Table 5. Inferred intrinsic parameters for shape (axis ratios q\ and 92) 
and orientation (cosine of the inclination angle cos Central values and 
dispersions are the mean and the standard deviation of the PDF. 



< ±8/nV s 2 t per cent. We consider a further error from the pecu- 
liar modelling we used. We looked for the minimum of Xx+sz by 
modelling only the data points beyond a central cut of 100 kpc with 
a simple /3-model for the density profile. The shift in the minimum 
was Ae A < 0.02. 

We added all the discussed independent sources of error in 
quadrature and ended up with an additional uncertainty of Se B ^ B = 
0.11. 



7 DEPROJECTION 

The combined X-ray plus SZ analysis allow us to infer the width 
of the cluster in the plane of the sky (parameterized in terms of the 
ellipticity e) and its size along the line of sight (expressed as the 
elongation eA). We have to use these two observational constraints 
to infer the intrinsic shape of the cluster (q± and 32) and its orien- 
tation ($ and 4>). The problem is clearly under-constrained dSerend 
l2007t) . The use of some a priori hypotheses on the cluster shape 
can help to disentangle the intrinsic degeneracy. We apply here 
some Bayesi an methods already employed in grav i tational lens- 
ing analyses dOguri et al 1 20051 : ICorless et alj|2009l ISereno et al.l 
l2010al : ISereno & Umets J201 ll) . 

As likelihood function we exploited 



£ = 



°<T t 



■ exp 



(27T) 1 /- 

P(eA — AeA) exp I — 



2<r? 
1 



(21) 



Ae A 



Ellipticity e and elongation ba are functions of gi, 52, 1?, and 
<f>. -P(ca) is the marginalized posterior probability distribution 
for the elongation parameter obtained in Sec. [6] The distribution 
was smoothed using a Gaussian kernel estimator dVio et al.|[l994l : 
Rvden 1996). We parameterized the additional statistical and sys- 
tematic uncertainty on the elongation as a shift following a nor - 
mal distribution with dispersion equal to 5e B ^ B dD'Agostinj||2003l) . 
Since the systematic error is quite smaller than the dispersion the 
impact on the final results is negligible. 

As prior for the intrinsic shape, we considered a flat distri- 
bution for the axial ratios in the range g m i n < q% ^ 1 and 
<?i ^ 92 ^ 1. Probabilities are defined such that the marginal- 
ized probability P(qi) and the conditional probability P(g2|gi) are 
constant. The probabilities can then be expressed as 



P(«i) = 1/(1 - 9mm) 

for the full range g m i n < gi ^ 1 and 

p(qz\qi) = (1 - Qi)' 1 



(22) 



(23) 



for g2 ^ qi and zero otherwise. The resulting probability for g2 is 
then 



P(<?2) = 



1 - gmin 



hi 



1 - gmin 

1 - qi 



(24) 
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Figure 8. PDFs for the intrinsic parameters under the prior hypotheses of flat q-distribufion and random orientation angles. Panels from the left to the right are 
for qi, q2, and cosi9, respectively. The posterior PDFs are plotted as white histograms. Thin lines in the q-panels denote the a priori flat (/-distribution. The 
thin and the thick line in the $ panel represent a priori random or biased orientation, respectively. 
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Figure 9. Contour plot of the marginalised PDF for the the axial ra- 
tios 51 and c/2- Contours are plotted at fraction values exp(— 2.3/2), 
exp(— 6.17/2), and exp(— 11.8/2) of the maximum, which denote con- 
fidence limit region of 1-, 2- and 3-cr in a maximum likelihood framework, 
respectively. Thick lines denote the loci of points compatible with either 
prolate or oblate configurations. 
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Figure 10. Contour plot of the marginalised PDF for the the axial ra- 
tio qi and the inclination cos Contours are plotted at fraction val- 
ues exp(— 2.3/2), exp(— 6.17/2), and exp(— 11.8/2) of the maximum, 
which denote confidence limit region of 1-, 2- and 3-cr in a maximum like- 
lihood framework, respectively. 



A flat distribution allows also for very triaxial clusters (c/i < qi <C 
1), which are preferentially excluded by iV-body simulations. We 
fixed c/ m i n = 0.1. 

For the orientation, we considered a population of randomly 
oriented clusters with 



p(cos$) = 1 

for ^ cos ■£) ^ 1 and 

P(0) = - 



(25) 



(26) 



for -tt/2 < <f> < 7r/2. 

The posterior distribution were investigated by running four 
Markov chains and checking for convergence. Results are summa- 
rized in Table [5] and Fig. [8] Inferred distributions are reassuringly 
dominated by the likelihood. The prior on the axial ratios plays an 
heavy role only for the final distribution of r/ 2 , which is not well 
constrained by data. The distribution of qi is peaked at ~ 0.8 with 



a tail in correspondence of more triaxial shapes. The spherical hy- 
pothesis is ruled out by the very precise measurement of the pro- 
jected ellipticity. Even if a priori orientations are random, a pos- 
teriori biased orientation are favoured. In Fig. [8] we compare the 
inferred distribution for the orientation of the major axis with re- 
sult s from semi-analytic al dOguri & Blandforl l200^) and numeri- 
cal ( Hennawi et al J2007I) investigations, which showed a large ten- 
dency for lensing clusters alike A 1689 to be aligned with the li ne 
of sight. Such condition can be expressed as dCorless et all 2009) 



p(cosi?) oc exp 



(cos 1? 



2*3 



(27) 



with a$ — 0.115. Noteworthily, the inferred distribution for cosi? 
is quite similar to what expected for biased clusters. 

Two dimensional posterior probabilities can provide further 
insights, see Figs. |9l and 1 101 Mildly triaxial shapes (0.7 < qi < 
t/2 < 0.9) are favoured, see Fig. [9] Oblate configurations (c/2 = 1) 
are generally excluded, apart from a small region near qi ~ 0.8. On 
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the other hand, prolate shapes (gi =92) can reproduce the obser- 
vational constraints. The major axis is preferentially aligned with 
the line of sight (cosi9 = 1), see Fig.QJJ] Orientations in the plane 
of the sky (cosi? = 0) or intermediate inclinations (cos # ~ 0.5) 
are compatible with mildly triaxial shapes (gi ~ 0.8). Very triaxial 
shapes, qi > 0.3, are associated only with aligned configurations 
(cosi?<l). 



8 CONCLUSIONS 

Knowledge of the intrinsic shape and orientation of halos is crucial 
to unbiased determinations of their masses and concentrations and 
to infer their hydrodynamical status. We have discussed a method 
to determine the geometrical properties of the ICM in rich galaxy 
clusters. The method is direct and the only hypothesis is that the 
cluster halo is approximately ellipsoidal. No assumption is needed 
about hydrostatic equilibrium. We exploited combined information 
from deep X-ray and SZe observations to constrain the elongation 
of the cluster gas distribution. This measurement of the size along 
the line of sight can be used with the measured width in the plane 
of the sky to infer the intrinsic form. Since the problem is under 
constrained, Bayesian inference has to be deployed. We applied 
the method to A1689 and found that the gas distribution is pref- 
erentially triaxial and elongated along the line of sight, in agree- 
ment wkli_corripJ^meirtaryrecent results derived with lensing anal- 
yses. ISereno & UmetsvJ d201ll) found a minor to major axis ratio 
~ 0.5 ± 0.2 f or the shape of the m atter distribution. Similar values 
were found in lCorless etai] d2009h and lMorandi et alj ( 1201 ll) . 

The technique described in the present paper provides direct 
estimates on shape and orientation by exploiting very simple as- 
sumptions. It does not rely on the hydrostatic equilibrium hypothe- 
sis and we could infer the 3D form without any assumption derived 
from iV-body si mulations. The me t hod is a de velopment of t he ide a 
first presented in |Pe Filippis etalld2005l) and lSereno et al.1 d2006l) . 
We implemented some significant improvements, i) The method is 
still parametric but does not rely anymore on the simple isothermal 
j3 model. The employed profiles can mimic complex features in ei- 
ther the electronic density or the temperature profile, ii) Instead of 
the central Compton parameter yo, we considered the more reliable 
integrated Compton parameter. Hi) Ev en if astronom ical deprojec- 
tion is an under-constrained problem JSerenoll2007t) . we could in- 
fer the 3D structure of the cluster without assuming any specific 
configuration. Our statistical analysis relied on Bayesian methods 
which enable to make an inference about a number of variables 
larger than the number of th e observed data. On the other hand in 
|Pe Filippis et alj d2005l) and lSereno etal.1 (|2006), the 3D distribu- 
tion was assumed to be either triaxial and aligned with the line of 
sight or prolate or oblate, iv) On top of this, we used more recent 
data which allowed us to infer the density and temperature profiles 
up to r < 10 3 kpc. 

Thanks to the XMM spectroscopic data and the Chandra de- 
rived surface brightness profile, X-ray and SZe analyses probe sim- 
ilar regions. The temperature profile obtained from Chandra is con- 
sistent in shape with the XMM-MOS results, but is systematically 
higher by about 1 1 per cent, on ave r age. S imilar results are reported 
in the study by iNevalainen et alj d2010t) of the cross-calibration 
uncertainties in the X-ray measurements of galaxy clusters. Fur- 
thermore, thanks to the larger radial extension of XMM observa- 
tions, we could detect a temperature decrement at larger radii. The 
smaller temperature and the decrement at large radii account for the 



gap between the estimate o f eA ~ 0-7 in the prese nt paper and the 
rounder value of eA ~ 1 in Pe Filip pis et al. I d2005h . 

The application of the method to A 1689 strengthens the view 
of a triaxial cluster elongated along the line of sight. The estimated 
triaxial structure of A 1689 together with the inferred orientation 
offers a coherent scenario explaining at the same time lensing, X- 
ray and SZ effect observations. A 1689 emerges as a very massive, 
slightly triaxial halo elongated along the line of sight, with a con- 
centration just above what predicted from numerical simulations 
and not so far from hydrostatic equilibrium. 

A triaxial model of the dark matter halos based on the a nalysis 
of iV-body simulation s was proposed bv lJing & Sutol (2002). Later 
on. lLee&Sutol ( l2003l) considered the relation between the gas dis- 
tribution and the matter profile in this model under th e hypothe- 
sis of hydrostatic equilibrium. The model was tested by Kawahara 
(2010), who found that the observed distribution of projected axis 
ratios of 70 clusters imaged by the XMM-Newton satellite was con- 
sistent with predictions. 

According to the results in ljing & Sutol d2002l) and lLee & Sutol 
d2003l) . the intracluster gas distribution in hydrostatic equilibrium 
under the gravity of triaxial dark matter halo with g™ attcr = 0.5 ± 
0.2, as inferred from lensing dSereno & Umets J201 lh . should have 
qi ~ 0.8 ± 0.1. As expected, the gas distribution is rounder than 
the total matter density and is in good agreement with our estimate 
of gi =0.70 ±0.15. 

The central mass discrepancy between lensing and X-ray esti- 
mates assuming hydrostatic equilibrium might be solved by the ori- 
entation bias. Effects of tr iaxiality might pla y a larger role than de- 
viations from equilibrium. iPeng et al.l d2009) showed that masses as 
estimated with lensing and X-ray analyses might agree if the matter 
halo of A 1689 was oriented along the line of sight with an elonga- 
tion of ~ 0.6. Thanks to our combined X-ray and SZ analysis, we 
derived a direct estimate of the elongation of the gas distribution 
of eA ~ 0.7 ± 0.2. Being the gas rounder than the mass distri- 
bution, our estimate of the elongation is fully consistent with the 
value required to reconcile the discrepant estimates. On the same 
page,[Morandi et"al] J20 1 lh relying on theoretical extrapolations of 
results from iV-body simulations guessed from their combined X- 
ray plus lensing analysis that the non thermal contribution to the 
total pressure in A 1689 was about 20 per cent. 

Finally, our results reduce the over-concentration problem for 
A1689. Lensing analyses performed under the hypothesis of spher- 
ical symmetry derived an unusually high concentration for the 
dark matter halo of A1689, well in excess of predictions from N- 
body simulation. Triaxiality was proposed to reconcile observa- 
tions and theoretical pre dictions (Ogurietal. 2005; Coiie sTet al.l 
2009; ISereno&Umetsull201ll) . The shape and orientation of the 
gas distribution of A 1689 rese mble those derived with lensing anal- 
yses of the dark matter halo ( Sereno & Umetsvj|20Tlh . In particu- 
lar, we retrieved and confirmed the orientation bias that very likely 
causes the over-concentration problem. The halo of A 1689 seems 
to be elongated along the line of sight and only slightly over- 
concentrated but still consistent with theoretical predictions. 
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